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Abstract 

We consider viscous rotating accretion flows in which most of the viscously dissipated 
energy is stored as entropy rather than being radiated. Such advection- dominated flows 
may occur when the optical depth is either very small or very large. We obtain a family 
of self-similar solutions where the temperature of the accreting gas is nearly virial and 
the flow is quasi-spherical. The gas rotates at much less than the Keplerian angular 
velocity; therefore, the central stars in such flows will cease to spin up long before they 
reach the break-up limit. Further, the Bernoulli parameter is positive, implying that 
advection-dominated flows are susceptible to producing outflows. Convection is likely in 
many of these flows and, if present, will tend to enhance the above effects. We suggest 
that advection-dominated accretion may provide an explanation for the slow spin rates of 
accreting stars and the widespread occurrence of outflows and jets in accreting systems. 

Subject headings: accretion: accretion disks — black hole physics — hydrodynamics 



1. Introduction 

In astrophysical accretion flows, gravitational potential energy is converted to kinetic 
and thermal energy of the accreting gas. If the thermal energy is efficiently radiated away, 
the orbiting gas becomes much cooler than the local virial temperature and takes up a 
thin disk-like configuration. Accretion disk solutions with these characteristics, derived 
originally by Shakura & Sunyaev (1973) and Lynden-Bell & Pringle (1974), have been 
used to model a variety of accreting systems in astrophysics (see Frank, King, & Raine 
1992 for a review). 

We ask here a simple question: What is the nature of the flow if the accreting gas is 
unable to cool efficiently? We use the term advection-dominated to refer to a flow where 
the bulk of the liberated thermal energy is carried in by the accreting gas as entropy 
rather than being radiated. In this Letter we write height-averaged equations to describe 
advection-dominated flows and discuss the properties of a particular class of self-similar 
solutions. 



2. Height- Averaged Equations of Advection-Dominated Accretion 

As in the standard theory of thin accretion disks, we vertically average the flow equa- 
tions and consider a two-dimensional flow in the equatorial R(p plane. We assume a steady 
axisymmetric flow so that d/dt = d/d<p = and all flow variables are functions only of R. 
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Assuming gas pressure to dominate, we write the pressure as P = pc 2 where p(R) 
is the height-averaged density and c 2 (R) is the isothermal sound speed. We denote the 
Keplerian angular velocity by Qk(R) = (GM/ R 3 ) 1 / 2 and define the Keplerian velocity 
vk = (GM I 'R) 1 / 2 , where M is the central mass. The surface density of the gas is E = 2pH, 
where H ~ Rc s /vk is the vertical scale height. Following Shakura & Sunyaev (1973) we 
take the kinematic coefficient of shear viscosity to be v = ac s H = ac 2 /Q,K, where a is 
a constant. Purely on dimensional grounds, this approximation seems to be particularly 
appropriate for the self-similar flows discussed here. 

The density of the gas p, its radial velocity v, angular velocity O, and isothermal sound 
speed c s , satisfy the following four differential equations, namely the continuity equation, 
the radial and azimuthal components of the momentum equation, and the energy equation 
(e.g. Abramowicz et al. 1988, Narayan & Popham 1993): 
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In equation (4), the left-hand side is the advected entropy, where T is the temperature and 
s the entropy, while the right-hand side gives the difference between the energy input per 
unit area due to viscous dissipation (Q + ) and the energy loss through radiative cooling 
(Q~). For convenience we have defined a parameter e = (5/3 — 7) /(j — 1), where 7 is 
the ratio of specific heats; note that e = in the limit 7 = 5/3 and e = 1 when 7 = 4/3. 
Substituting the viscous dissipation rate for Q+, we obtain 

+ ^_ _ 2apc 2 s R 2 H fdtt\ 2 ^_ _ t 2apc 2 s R 2 H ( dQ\ 



The parameter / measures the degree to which the flow is advection-dominated. In the 
extreme limit of no radiative cooling, we have / = 1, while in the opposite limit of very 
efficient cooling, / = 0. Finally, we define e = e/f. This is an important parameter which 
plays a critical role in determining the nature of the flow. 

Let us for simplicity assume that e is independent of R. Equations (l)-(4) then 
permit a self-similar solution of the form 
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The density p may be obtained from the mass accretion rate, M = —AirRHvp. In equations 
(7)-(9), the first relation gives the exact solution, while the second corresponds to the limit 
when « 2 <C 1. If we include bulk viscosity through an a prescription there are corrections 
of order a 2 in the results. Also, if e < 0, i.e. if 7 > 5/3, there is a second class of self-similar 
rotating wind solutions where v > 0. We do not discuss these extensions here. 



3. Properties of the Self-Similar Solution 

In the limit of very efficient cooling, / — >• and e — >• 00, and the solution given 
in equations (7)-(9) corresponds to a standard thin accretion disk with v, c s <C vk and 
O — > VLk- The properties of such disks have been widely discussed in the literature. In 
this Letter we are interested in the opposite limit of advection-dominated flows where / is 
a reasonable fraction of unity and e ~ e < 1. In this limit, the self-similar solution derived 
above has several interesting properties which we now discuss. Although our discussion is 
based on this particular solution, we believe that the results are valid for general advection- 
dominated accretion flows. 

Equation (9) shows that the sound speed in advection-dominated flows is comparable 
to the Keplerian speed vk, which means that the temperature of the accreting gas is nearly 
virial. This is of course expected since by assumption the gas has no way to cool. When 
c s ~ vk, the disk vertical thickness H is comparable to R and the flow is quasi-spherical. 
Unfortunately, the vertical averaging on which equations (l)-(4) are based is suspect when 
H/R ~ 1, and therefore some of the numerical factors in our results may be in error. 

The radial velocity of the accreting gas is proportional to a. This is because the radial 
speed is determined principally by how fast the viscosity can move angular momentum 
outwards. Since v ~ ac 2 /vx, the radial velocity tends to be much larger in advection- 
dominated flows than in thin disks. 

When e is small, we see that the angular velocity Q of the flow is smaller than the 
local Keplerian by a factor ~ (e ) 1 / 2 . Indeed, as e — > 0, which corresponds to 7 — > 5/3, 
Q, goes to zero and our solution matches on to the Bondi spherical accretion solution for 
this 7. (In a sense, the solution (7)-(9) is a natural extension of self-similar 7 = 5/3 Bondi 
accretion to rotating flows and general 7.) Because advection-dominated accretion has a 
significantly sub-Keplerian O, if such a flow were to be present around an accreting star, 
the star would cease spinning up long before the rotation rate approached the break-up 
limit. This may have important implications for the spin rate of accreting stars. 

Equation (2) shows that the self-similar solution satisfies the relation 
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From this we can compute the normalized parameter b = Bejv 2 K , where Be is the Bernoulli 
constant: 

, ( 7 5\c^ = 3e-e^ 

2v 2 K + U-l V A 5 + 2e " 

As is well-known, Be is conserved in adiabatic flows in the absence of viscosity. Therefore, 
whenever b is positive, it implies that if we were somehow to turn the gas around and to 
let it flow adiabatically outward on a radial trajectory, the gas would reach infinity with a 
net positive energy. On the other hand, if b < 0, the gas cannot spontaneously escape to 
infinity. Equation (12) shows that b is positive in advection-dominated flows for all values 
of a and for any 7 < 5/3 so long as / > 1/3. Winds, jets and other outflows are common 
in several accreting systems in astrophysics and it is tempting to speculate that these may 
originate in advection-dominated flows. 

Note that the positivity of the Bernoulli constant does not imply a lack of conservation 
of energy, but simply arises because viscous stresses transfer energy from small to large 
radii. Even in the standard theory of thin accretion disks, the energy radiated from any 
given annulus of the disk is greater by a factor of 3 than the net gravitational energy 
released within that annulus (e.g. Frank et al. 1992). A similar phenomenon occurs in our 
solutions, and the excess b at any radius represents energy transferred to that radius from 
smaller radii. When proper boundary conditions are applied and the full global problem 
is solved, total energy will of course be conserved. 

It is easily shown that the entropy increases inwards in our self-similar solution. We 
need therefore to consider the possibility of convective energy transport. In a rotating 
medium, the condition for a dynamical convective instability is 

N *ff = N +K =-pdR-dR +K <0 ' (13) 

where N is the usual Brunt- Vaisala frequency and k is the epicyclic frequency which in 
our case is equal to O. For our self-similar solution we find (in the limit when a 2 <C 1) 

N 2 eff = )° e ' +6e /~ 15e ^K- (14) 
e// (5 + 3e)(5 + 2e') K V ; 



This gives a dynamical instability for 
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Even when / is below this limit a mild double diffusive instability will be present. 

Although advection-dominated flows are convectively unstable, we can show that con- 
vection, if at all, only strengthens our main conclusions. The height-integrated convective 
energy flux F c in the presence of an unstable entropy gradient can be written quite gener- 
ally as 

F c = ~ZK C T^, (16) 
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where K c is the diffusion constant associated with convective transport. In the spirit of 
the Shakura-Sunyaev viscosity formula let us write K c = a c c 2 s /VtK- Substituting our self- 
similar solution into equation (16) we then find that F c scales as R~ 2 and has a negative 
divergence. Convection therefore provides an additional energy source in equation (4) of 
magnitude 

tt = -*-#=-^T*. (17) 

This term has the same structure as the left-hand side of equation (4). We can therefore 
combine the two terms by simply redefining e to be 

e /- a c c 2 s \ e ( 2 a c 
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With this definition, all of the results in equations (6)-(10) continue to hold. 

In general we expect a > a c . This is because convective turbulence contributes to the 
viscosity but not all sources of viscosity produce bulk energy transport. (Ryu & Goodman 
(1992) show that if the entropy in an accretion disk is stratified vertically, then convec- 
tion actually moves angular momentum inwards; in our case we have a radially stratified 
medium and so we expect convection to produce a normal outwardly-oriented angular mo- 
mentum flux.) Since a > a c , equation (18) shows that convection only introduces a modest 
perturbation to the energy equation (4). This is very different from the situation in stellar 
interiors where convection, once it sets in, strongly modifies the structure and essentially 
forces an isentropic radial distribution. In our case, the strong advection ensures that con- 
vection can only be a moderate perturbation. Indeed, as equation (18) shows, the primary 

effect of convection is to reduce the value of the parameter e . This means that the various 
properties of advection-dominated flows discussed earlier become even more pronounced. 
In particular, when there is convection, the parameter b remains positive for a wider range 
of / than in the non-convective case. The reason is that convection transports energy from 
smaller to larger R and therefore adds to the effect of the viscous stress. 

The above discussion suggests that the main results of this paper are likely to survive 
even if there are other instabilities in the flow. 



4. Discussion 

Advection-dominated accretion is expected to occur in several astrophysical situations. 
Thin accretion disks for example exhibit a thermal instability at sufficiently low optical 
depth, when the cooling through free-free emission is unable to keep up with the viscous 
energy generation. Numerical models of accretion disks in cataclysmic variables (Narayan 
& Popham 1993) show that, under the influence of this instability, the disk switches to an 
advection-dominated flow where only a fraction of the released energy is radiated. A more 
extreme possibility is that at very low accretion rates the infalling material may never cool 
sufficiently to collapse to a thin disk, and we could imagine an advection-dominated flow 
all the way from the outermost radius down to the central star or black hole. 

Advection can dominate also in the opposite limit when the mass accretion rate is 
very high and the optical depth r of the disk is large. If the cooling time of the disk, 
t coo i ~ Ht/c, is longer than the accretion time R/v, most of the accretion energy is 
retained in the gas and we have an advection-dominated flow. Such conditions are possible 
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in young stellar objects and symbiotic stars in outburst (Popham et al. 1993) and in 
accreting neutron stars and black holes at high mass accretion rates. 

Thus, advection-dominated accretion is likely to occur in a number of astrophysical 
systems, and some objects may have both advection-dominated and standard cooling- 
dominated zones at different radii. Surprisingly, very few studies in the literature have 
included advection effects, and the physics of advection-dominated accretion has been 
hardly discussed at all. The self-similar solutions we present in §2 of this Letter represent 
a first step towards understanding such flows. 

Our self-similar solutions have some interesting features which we have described in 
§3. Here we would like to highlight the following: 

1. The angular velocity O of the flow is less than the Keplerian angular velocity VL K ; in 
fact, for the astrophysically interesting case of 7 — > 5/3 (e — > 0), we find O <C Ok- In 
accretion theory it is usually assumed that the infalling material rotates with O ~ 
and that accreting stars will spin up to the "break-up limit" corresponding to a stellar 
spin rate fi* ~ Qk(R*), where R* is the equatorial radius of the star. Advection- 
dominated accretion flows are very different and the star may reach an equilibrium 
spin state with fi* <C Qk(R*)- At the equilibrium fi* the system will adjust the 
angular momentum flux in such a way that remains constant (Popham & Narayan 
1991, Paczyhski 1991). Similarly, in the case of accretion onto a magnetic star it 
is usually assumed that the central star spins up until fi* = Qk(Ra), where Ra is 
the Alfven radius. Once again, if the accretion is advection-dominated, we can have 
fi* < Q K (R A ). 

2. The scaled Bernoulli parameter b (see equation 12) is positive in self-similar advection- 
dominated flows for a wide range of parameters. We have argued that this means the 
gas in such flows is capable of spontaneously escaping to infinity. We suggest that 
advection-dominated flows may provide a generic explanation for many of the outflows 
and jets that are so ubiquitous in accreting systems. It is beyond the scope of this 
Letter to discuss the exact mechanism by which outflows may be generated. One 
possibility is that a shock may divert part of the high entropy accreting material 
into an outward-pointing trajectory. Alternatively, the positivity of b may make the 
generation of winds so easy that even a quite modest radiative or magnetic stress may 
set off a substantial wind. (We must caution however that, in some circumstances, 
once the material reaches the surface of the disk it may cool and lose its positive 6, 
perhaps making it harder to generate a wind.) 

3. By definition, the radiative luminosity of an accretion-dominated flow is much less 
than the standard GMM/R* that is usually associated with accretion. Therefore, 
estimates of the mass accretion rates of observed systems could be seriously in error, 
especially when the central star is a black hole which can swallow the advected entropy, 
or when much of the accretion energy escapes as an outflow. Moreover, since the 
accreting gas is almost at virial temperature, the spectrum is likely to be much hotter 
than the spectrum of an equivalent thin disk. This is particularly true of optically 
thin advection-dominated flows. 

4. Advection-dominated flows are convectively unstable. Because convection transfers 
energy from small to large radii, we find that it enhances the effects described above. 
The self-similar solutions described here technically extend from R = to R = 00 with 

no boundaries. How relevant are the results if the accretion flow is advection-dominated 
only over a finite range of R and is bounded on the outside and inside by non-advection 
zones? We have carried out some preliminary investigations of this question by numeri- 
cally solving equations (l)-(4) for various boundary conditions. Figure 1 shows an example 
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where the accreting gas starts off as a thin disk at large R, goes through an advection- 
dominated phase, and becomes a rotating settling "star" at small R. Over a range of 
intermediate radii the numerical solution is close to the self-similar form. We have calcu- 
lated other numerical solutions where the central star is a black hole and the flow goes 
through a sonic point, and again we find that the flow is nearly self-similar in between the 
two boundaries. 

From these numerical examples we conclude that the self-similar solutions we have 
described may be of more than academic interest and that perhaps real flows often re- 
semble the self-similar form. In any case, the main features of these solutions, namely 
slow rotation and the ability to make outflows, are probably generic to a wide class of 
advection-dominated flows. 
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Figure Caption 

Figure 1. A numerical solution of equations (l)-(4) with GM = 1, a = 0.1, / = 1, 
and e = 1/3 (or 7 = 1.5). At R = 10 6 (the outer edge), the boundary conditions are 
O = Qk and c s = 0.001i>k, appropriate to a very thin disk. At R = 1 (the inner edge), 
the boundary condition is v = 0.05c s , corresponding to a settling star. The "star" in this 
example extends out to about log(R) ~ 1. The angular momentum flux has been set to 
zero, corresponding approximately to an equilibrium spinup state. (The precise equilibrium 
value of this parameter depends on the structure of the central star, cf Popham & Narayan 
1991). The solid lines in the three panels show the radial velocity v, the scaled angular 
velocity VL/VLk, and the scaled Bernoulli parameter b. The dashed lines show the results 
corresponding to the self-similar solution. 
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